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Abstract 

An aspect of fluid dynamics lies in the search of possible statistical models for Navier-Stokes (NS) 
fluids described by classical solutions of the incompressible Navier-Stokes equations (INSE). This 
refers in particular to statistical models based on the so-called inverse kinetic theory (IKT) . This 
approach allows the description of fluid systems by means a suitable 1-point velocity probability 
density function (PDF) which determines, in terms of suitable "moments", the complete set of fluid 
fields which define the fluid state. A fundamental related issue lies in the problem of the unique 
construction of the initial PDF. The goal of this paper is to propose a solution holding for NS 
fluids. Our claim is that the initial PDF can be uniquely determined by imposing a suitable set of 
physical realizahility constraints. 

PACS numbers: 05.20Jj,05.20.Dd,05.70.-a 



I. INTRODUCTION: CSM-INSPIRED STATISTICAL MODELS 



Fundamental aspects of fluid dynamics are related to construction of statistical models 
{ /, r} for fluid systems. These are sets {/, F} formed by a suitable probability density 
function (PDF) and a phase-space F (subset of M") on which / is deflned. By deflnition, a 
statistical model {/, F} of this type must permit the representation, via a suitable mapping 

{f,r}^{z}, (1) 

of the fluid fields {Z} = {Zi, i = l,n} which deflne the state of the same fluids. Depending 
whether the mapping provides the complete set or only a subset of {Z} the statistical model 
is denoted respectively complete or incomplete. In the case of a Navier-Stokes (NS) fluid, i.e., 
an incompressible isothermal and isentropic Newtonian fluid described by the incompressible 
NS equations (INSE), the fluid flelds are 

{Z}^{Po,\,p,,St}, (2) 

where in particular po (the mass density) and St (the thermodynamic entropy) are both 
assumed constant in x / (with denoting the closure of the conflguration domain i7 C 
and / C M a time interval/axis). Moreover, V and pi denote the fluid velocity and the 
(strictly positive) kinetic pressure 

p,{r,t)=p{r,t)+po{t) + <p{r,t), (3) 

where p{r, t),pQ{t) and 0(r, t) represent respectively the fluid pressure, the (strictly-positive) 
pseudo-pressure and the (possible) potential associated to the conservative volume force 
density acting on the fluid. The construction of a statistical model involves the introduction 
of a suitable PDF of some sort (/), deflned on an appropriate phase-space F in terms of which 
the fluid flelds (or only a proper subset of them) can be represented via suitable velocity or 
phase-space moments of the PDF. A well-known example of incomplete statistical model. 



holding for incompressible NS 
developed originally by Hopf 



iuids, is provided by the so-called statistical hydromechanics 
it and later extended by Rosen and Edwards {HRE 
approach), based on the statistical model {fH,^i}- This relies on the introduction of the 
1— point (or local) velocity- space PDF, /i, to be intended as the conditional PDF of the 
velocity v {kinetic velocity) with respect the remaining variables. In the HRE approach these 
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are identified with (r,t), where (r,t) E Q x I, while /i = /i(r, u,t; Z), with u = v — V(r,t) 
the relative kinetic velocity, is identified with 

/^^5(v-V(r,t)) (4) 

fn denoting the three-dimensional Dirac delta defined in the velocity space f/ C M'^, with v 
belonging to U and (r,t) E Qx I (with C and I the configuration domain and the time 
axis). Hence, fn is defined by assumption on the set spanned by the state vector x =(r, v), 
i.e., the restricted phase space Ti = Q x U . It follows that, in this case, only one velocity 
moment of /i, corresponding to G = v, is actually prescribed in terms of the fluid fields and 
reads 

j rf3vG/H(r,u,t;Z)=V(r,t). (5) 

The statistical model {/h, Fi}, which was adopted also by to Monin {4] and Lundgren j^, 
belongs actually to a more general class of statistical models inspired by Classical Statistical 
Mechanics (CSM). In this context, a convenient alternative approach - which permits instead 
the representation in terms of / of the complete set of fiuid fields (hence providing a com- 
plete statistical model) - is provided by with the so-called IK T- statistical model {/i,ri}, 
developed in the framework of inverse kinetic theory (IKT; Tessarotto et al. 7Hll|). This is 
based on the construction of a suitable 1-point velocity-space PDF, /i, defined in such a way 
to yield,, in terms of a prescribed set of velocity moments of the same PDF, the complete 
set of fluid fields characterizing a NS fiuids. 

A basic issue is therefore represented by the construction of the general solution fi (t) for 
the IKT-statistical model, and in particular the initial PDF fi{to), based on the available 
information (i.e., observables) on the fiuid system. 

The goal of this paper is to point out that /i(t) can be uniquely determined at the initial 
time t = to- The basic results are represented by THMs 1 and 2 (see Sections 4 and 5). 
We claim that, subject to the requirement of factorization at the initial time for all multi- 
point PDFs the initial PDF /i(to) can be determined imposing the principle of entropy 



12 1 . In detail, we intend to 



maximization (PEM [13|) on the Boltzmann- Shannon entropy 
show that the initial PDF /i(to) can be determined in such a way to satisfy the following 

requirements: 

1. (requirement #1) it satisfies the complete set of constraints placed at the initial time 
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t = to by the observables defined for the fluid system (physical realizahility conditions)'^ 

2. (requirement #2) it does not require the specification of liiglier-order PDFs, witli 
s > 2; 

3. (requirement #3) it provides the general solution for the initial PDF required by the 
IKT- statistical model {/i, Fi} . 

II. PHYSICAL OBSERVABLES OF A NS FLUIDS 

An important preliminary task to accomplisli is to determine tlie complete set of physical 
observables required to define the functional class of the PDF {{fi(t))^} , to be prescribed 

le complete set of the fluid field which 



ll| for an isothermal, isentropic and 



at the initial time to- These include, in particular, t 
specify the state of the fluid. In the IKT approach [7 
incompressible NS fluid the latter must be represented in terms of the set 

A further observable can, however, be defined in terms of the velocity field V (r,t) . This is 
represented by the continuous velocity- frequency density (C-VFD), f['^^^'^\t,v), namely the 
velocity-frequency associated to the fluid velocity occurring in subsets of the fluid domain 

C M^, to be defined in particular so that 

d\ff''^'\t,^r) = 1, (6) 
u 

U = M.^ denoting the three-dimensional velocity space. For definiteness, let us assume in 
the remainder that: 

• Assumption #1: the fluid domain fi is a bounded subset of with finite measure 
(for example is a cube). 

• Assumption #2: the fluid velocity V(r,t) is bounded. Hence there exists a positive 
constant Vb G M, with VB(t) = 2 sup {Vfc(r,t); = 1, 2, 3; r G fi} and we can define the 
bounded subset of 



f/ = <; V : \vA < -Vr; V gM^ 



to be identified with the velocity space of the statistical model. 
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/r^^(t,v)= lim /ir(t,v), (7) 



Let us require for simplicity that, thanks to Assumption #1, Q is partitioned in hke 
cubic cells. Then let us introduce the notion of continuous velocity- frequency density (C- 
VFD) fi'^^^''\t,v) to be defined as the limit function 

/^■^^''■'(t, v) indicating the corresponding (discrete^ velocity- frequency density function (V- 
FDF) defined on a N-cell partition of Q 

/}S^''Ht,v) = l 5^iVi(r„v,t). (8) 

1=1, N 

To define A'^i(rj, v,t) let us notice that - thanks to Assumption #2 - for each component of 
the fluid velocity Vk{r,t) (with k = 1,2,3) the inequality 

\V,{r,t)\<lvB (9) 

holds. Hence, the velocity space U can be partitioned in M like cubic cells, with M G No 
denoting an arbitrary integer. Therefore, if denotes the position of the center of mass for 
the i — th configuration-space cell, A''i(ri, v,t) can be defined as the frequency of occurrence of 
the velocity fluid field v, assumed to belong to the velocity-space cell defined (for k = 1,2, 3) 
by the inequalities \Vk{ri,t) — Vk\ < j^- Thus, A''i(rj, v,t) can be defined as 

iVi(r„v,t) = ^ II e,fc(v), (10) 

fc=l,2,3 

e,,(v)^0(\4(r„t)-i;,-^) 

eivk-V,ir„t) + ^), (11) 
with B(x) denoting the Heaviside theta function, c G M the normalization constant 

^= / d'vj^ E n (12) 

^ i=l,Af fc=l,2,3 

and M, G Nq arbitrary integers. We stress that they can always be defined so that 
M = M{N) is a strictly monotonic function of (in particular, M and can be, for 
example, so that A^ = M^). Then, thanks to positions f|T0|) - f|T2|) . by construction fl^^^'^\\-,t) 
satisfies the normalization condition 

\\W;,^'\^r,t) = l. (13) 



This manifestly implies for the limit function f^-''^'^'^' ( v,t) the analogous condition of normal- 
ization (|6]). 

Additional observables, to be defined in analogy to f[-^^^'^\v ,t) , are however represented 



1. the 2-pomt velocity difference-frequency density function (VD-FDF) ^'^(r, v,to) to be 
identified with the frequency of the velocity difference v = Vi — V2 occurring between 
two positions ri and r2 with displacement r = ri — r2 and subject to the normalization 



u 

2. as well as the analogous s-point velocity-difference FDFs which can be defined in 
principle for arbitrary s > 2. 

It is obvious that in principle all the constraints provided by the infinite set of velocity- 
difference PDFs should he satisfied by the initial 1-point /i(to) [which defines the IKT- 
statistical model {/i,ri}]! 

III. PHYSICAL RE ALIZ ABILITY CONDITIONS FOR fi{to) 

In this reference, a natural question arises, i.e., whether the arbitrariness [in the defi- 
nition of {/i,ri}] can be used, hy proper prescription on its functional class {{fiito))^} , 
to determine it uniquely consistent not only with INSE hut also with the relevant physical 
ohservahles. An important preliminary task to accomplish is to establish the relationship 
of /i with the fluid fields. More precisely, here we state that, hesides the complete set of 
fluid fields evaluated at the initial time {V (r,to) ,Pi (i",io) , ST{to)} , the PDF must also he 
suitahly related to the initial 1— point velocity-frequency density function J\^^^'^\\,to)- 

In the following we shall require that fi satisfies the following constraints (to be intended 
as physical realizahility conditions for /i): 

Realizability condition ^1: it admits for all (r,t) G ^7 x / (including the initial time to) 
and G = 1, v,u^/2, uu, un^/2, Zn/i(r, v,t) the velocity and phase-space moments J dwGfi 



by: 




(14) 



u 
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and J dvfi In /i and satisfies the constraint equations (denoted as correspondence principle): 

J dvG/i(r,v,t) = l,V(r,t),pi(r,t), (15) 
Sim) = St, (16) 



u 



with S{fi{t)) = — J (ix/i(r, v,t) In /i(r, v,t) denoting the Bohzmann-Shannon statistical 
entropy associate to /i(t) and Fi the phase space Ti = Q x U; 

Realizability condition #2: at the initial time t = to it satisfies the constraint: 

{Mto))n = ff'''\^,Q- (17) 
where fi(t) = /i(r, v,t; Z) and {fi(t))^^ denotes the f2— average at time t 

n 

Realizability condition #3: at the initial time t = to, we require that the 2-point PDF 
/2(l,2,to) satisfies the constraint: 

Ur,^r,to) = ft\r,^r,to), (19) 
where /2(r, v,to) denotes the velocity-difference 2-point PDF defined as 

/2(r,v,t,) = -L.jdRj dVUl,2,to). (20) 

Then requiring that f2{l,2,to) is factorized 

Mi,2,to) = Mi,to)W,to) (21) 

this implies the constraint 

-L.JdRj dV/i(l,t„)/i(2,t„) = ^^'^^''(r, v^). (22) 

IV. CONSEQUENCES OF PEM 

Let us now show how fi{htri) = /i(xi, to; Z) can be determined for the IKT-statistical 
model {/i,Fi} by imposing |8|, |9| that at t = to it satisfies the constrained maximal vari- 
ational principle (also known as principle of entropy maximization or PEM; Jaynes, 1957 



13|: 

55(/i(l,t)) = 0, (23) 

together with the physical realizability conditions defined above. For definiteness, let us 
assume that /i(l,to) is an ordinary, strictly positive function of the general form 

with 

(/i(l,to))n 

to))j-^ determined by the constraint Eq.( !T7|) and h(to) to be assumed as a strictly 
positive and regular real function. 

Let us determine the initial PDF f to) subject to the whole set of constraints placed 
by the aforementioned by the physical realizability conditions. Invoking the factorization 
condition (12 ip and the position f l24p the extremal 1-point PDF /i(l,to) can be determined 
in an equivalent way either from the variational principle fl2^ of for the 2-point entropy 
principle 

5SiMto))=0. (26) 

In fact thanks to it follows that ^(/z) = -2/i(fi) / dxifi{l,to) ln/i(l,to) = 2/i(fi)5(/i). 
Hence one obtains 

SS{f2) = -2fi{Q) j dxMl,to)5h{l,to)[l + \nh{l,to)+ (27) 

+ Ao(l,to) + A2(l,to)M(l,to)' + j dX2\3{r,Y,to)fl{2,to) 







where A3(r, v,to) is a suitable Lagrange multiplier to be determined in such a way to ful- 
fill the constraint (12T]) . Then, considering variations of h{l,to) of the form 6fi{l,to) = 
fi{l,to)Sh{l,to), i.e., defined so that there results identically 5/i(l,to) = 0, the previous 
variational principle requires that h{l,to) must fulfill the variational equation 

/ d^,6h{l,to) yi^'^??" [1 + ln/^(l,t,) + Xoivuto) (28) 



+X2iri, to)u{l, to)'^ + j rfx2A3(r, v,to)/i(r2, V2,to) 



0, 
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with summation understood on the index i (for i = 1, 2, 3). The result, which provides the 
sought solution for /i(to), has the flavor of: 



THM.l - General solution for fi{to) 

The general solution of the PEM variational principle [Eg. ( f^) / in the functional class 
{(/i('^o))n} determined by realizability conditions #l-#3 is of the form [24\ ) with h{l,to) 
taking the form of a non- Gaussian distribution 

h{l,to) = exp |-1 - Ao(ri,t„) - A2(ri,to)M(l,to)2 - A3(l,t,)} , (29) 

where Xo{r,to) and A2(r,to) suitable Lagrange multipliers to be determined imposing at t = to 
the moment equations 17^) . while A3(l,to) denotes 

A3(l, to) = X-ii^i, vi, to) = j rfx2A3(r, v,to)/i(r2, V2, to) (30) 

and A3(r, v,to) is determined in such a way to fulfill the constraint fl22|) . 
It follows that: 

1. fi{l,to) is strictly positive and determined as velocity and phase-space moments the 
fluid velocity \_(ri,to), the kinetic pressure pi{ri,to) as well as the thermodynamic 
entropy St |8|410|; 



2. due to the arbitrariness of (/i(l,to))Q, /i(l,^o) is generally non-Gaussian and non- 
isotropic in velocity space; 

3. the PDF fi{to) does not require the specification of constraints on any higher-order 
multi-point PDF of order s > 2. 

PROOF 

In fact there it results identically 

6 [ rfxMt,) ^/^^f^''^°^,^^ ln ^£^,^^'^°^,^" =0. (31) 
J {h{l,to))^ {h{l,to))n ^ ' 

From the Euler-Lagrange equation it is immediate to reach Eq.f l29p for h{l,to) [see the 
Euler-Lagrange fl25]) above]. To prove 1) and 2) we notice that, by construction f to) 
is strictly positive its moments satisfy the correspondence principle constraints, while from 
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(12^ and (1291) it follows that /i(l,to) is generally non-Maxwellian. Finally to prove 3) it is 
sufficient to notice that constraints of the type f l22|) to be placed s-point PDFs (with s > 3) 
leave unchanged the Euler-Lagrange equation (1251) . Q.E.D. 

As a consequence of THM.l, Eqs. and fl2U]) provide the general solution for the 

initial PDF f\{l,to), fulfilling all the physical realizability conditions defined above [see 
Sec. 3]. Remarkably, the solution determined in this way does not require the specification 
of possible infinite set of additional constraints, analogous to Eq. (l22l) . to be placed on the 
initial multi-multipoint PDFs of order s > 2. 

V. PARTICULAR SOLUTIONS 

It is immediate to show that particular solutions are respectively provided by: 

a) the generally non-Gaussian PDF: 

with h{l,to) the Gaussian PDF 

h{l, to) = exp {-1 - \o{rio, to) - A2(rio, to)u{\, t^)^} ; (33) 

b) the Gaussian PDF 

f,{l,to) = /M(u(ri^);pi(ri,to)), (34) 

with 

fM{u{r,to);pi{r,to)) = , 3^ — rrexp < _ ^J^^^'^"^ i , (35) 

and Vthpiri,to) = \/2pi{ri,to, a)/po denoting the thermal velocity associated to the kinetic 
pressure. 

THM.2 - Particular solutions for fi{to) 

fi(to) coincides either with [3^) or ^^3^ respectively if : 

A) there results identically: 

A3(ri,Vi,to) =0; (36) 

B) besides ( f^) . also the equation 

= 1 (37) 

{h{to))n ^ ' 
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holds identically. 
PROOF 

In validity of the constraint (EH]) Eqs. ^ and dSS]) follow from Eqs.(I2S]) and ([30]) • 
In validity of Eq.f p7|) too, it is immediate to prove that PEM implies that for arbitrary 
variations 6h{r, u,to) it must result identically: 

J dx5h{r, v,to) {1 + In h{v, v,to)+ (38) 
+A„+A2u2} = 0. 

It follows that the Lagrange multipliers Ao(r,to) and A2(r,to), determined by imposing the 
correspondence principle, necessarily require 

fiito) = /a/(v - V{r,Q;pi{r,Q). (39) 

Q.E.D. 

VI. CONCLUSIONS 

An axiomatic approach, based on the IKT-statistical model {/i, Fi}, has been developed 
to determine the initial condition for the 1-point PDF fi which characterizes an isothermal 
and isentropic incompressible NS fluid. 

Lq_ particular, we have proven that - extending the statistical approach earlier developed 



11 j - the initial 1-point PDF can be uniquely determined by invoking the principle of 
entropy maximization. The present theory is based on the assumption that the initial PDF 
satisfles suitable physical realizability conditions. In detail, this requires that the initial 
PDF /i(to) be constructed in such a way to satisfy the physical requirements placed by : 

• the correspondence principle [in particular Eqs. ffTSj) ]. prescribing that the PDF deter- 
mines in terms of suitable moments the complete set of fluid flelds; 

• the velocity- frequency density function (V-FDF) J[^^^'^\\^to) f ll7p : 

• the 2-point velocity difference-frequency density function (VD-FDF) //'"'^''(r, v,)f:o), in 
term of the constraint ( 122|) . 

The theory has several important consequences: 
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1. THM.l provides the general form of the initial 1-point PDF satisfying the previous 
physical requirements; 

2. the initial 1-point PDF is generally non-Gaussian and non-isotropic in velocity-space 
(see THM. 1). 

3. particular solutions include, including the Gaussian 1-point PDF, are discussed in 
THM.2. 

In addition, the present theory provides an answer to the requirements placed in Sec.l, 
namely that the initial PDF /i(to) obtained in this way is unique. In fact we have proven 
(THM.l) that the determination of /i(to) does not require the specification of higher-order 
PDFs of order s > 2 and thus provides the general solution for the initial fi{to) of the 
IKT- statistical model {/i, Fi} . 
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